SDSS J123813.73-033933.0, a cataclysmic variable evolved beyond the period 

minimum. 

A. Aviles 1 , S. Zharikov 1 , G. Tovmassian 1 , R. Michel 1 , M. Tapia 1 , 

Institute) de Aatronomia, Universidad Nacional Autonoma de Mexico, Apartado Postal 877, 

22800, Ensenada, BC, Mexico 

M. Roth 2 , 

Las Campanas Observatory, Carnegie Instututio of Washington, Casilla 601, La Serena, Chile 

V. Neustroev 3 , 

Centre for Astronomy, National University of Ireland, Galway, Newcastle Rd., Galway, Ireland 

C. Zurita 4 , 

Instituto de Astrofisica de Canarias, c/ via Lactea s/n, La Laguna, E38200, Tenerife, Spain 

M. Andreev 5 , A. Sergeev 5 , 
Institute of Astronomy, Russian Academy of Sciences, Terskol, Russia 

E. Pavlenko 6 , 

Crimean Astrophysical Observatory, Nauchny, Ukraine 
V. Tsymbal 7 , 

Tavrian National University, Dep. Astronomy, Simferopol, Ukraine 

G.C. Anupama 8 , U.S. Kamath 8 , D.K. Sahu 8 
Indian Institute of Astrophysics - CREST, Bangalore 560 034, India 

1 Instituto de Aatronomia, Universidad Nacional Autonoma de Mexico, Apartado Postal 877, 22800, Ensenada, 
BC, Mexico 

2 Las Campanas Observatory, Carnegie Instututio of Washington, Casilla 601, La Serena, Chile 

3 Centre for Astronomy, National University of Ireland, Galway, Newcastle Rd., Galway, Ireland 

4 Instituto de Astrofisica de Canarias, c/ via Lactea s/n, La Laguna, E38200, Tenerife, Spain 

institute of Astronomy, Russian Academy of Sciences, Terskol, Russia 

6 Crimean Astrophysical Observatory, Nauchny, Ukraine 

7 Tavrian National University, Dep. Astronomy, Simferopol, Ukraine 

8 Indian Institute of Astrophysics, II Block Koramangala, Bangalore 560034, India 



- 2 - 



ABSTRACT 

We present infrared JHK photometry of the cataclysmic variable SDSS J123813.73- 
033933.0 and analyze it along with optical spectroscopy, demonstrating that the binary 
system is most probably comprised of a massive white dwarf with T e g- = 12000 ± 1000 K 
and a brown dwarf of spectral type L4. The inferred system parameters suggest that 
this system may have evolved beyond the orbital period minimum and is a bounce-back 
system. 

SDSS J1 23813.73-03393 3 .0 sta nds out among CVs by exhibiting the cyclical vari- 



ability that IZharikov et al.l ((2006) called brightening s. These are not related to specific 
orbital phases of the binary system and are fainter than dwarf novae outbursts, that 
usually occur on longer timescales. This phenomenon has not been observed extensively 
and, thus, is poorly understood. The new time-resolved, multi-longitude photometric 
observations of SDSS J123813. 73-033933.0 allowed us to observe two consecutive bright- 
enings and to determine their recurrence time. The period analysis of all observed 
brightenings during 2007 suggests a typical timescale that is close to a period of ~ 9.3 
hours. However, the brightenings modulation is not strictly periodic, possibly maintain- 
ing coherence only on timescales of several weeks. The characteristic variability with 
double orbital frequency that clearly shows up during brightenings is also analyzed. 

The Doppler mapping of the system shows the permanent presence of a spiral arm 
pattern in the accretion disk. A simple model is presented to demonstrate that spiral 
arms in the velocity map appear at the location and phase corresponding to the 2:1 
resonance radius and constitute themselves as a double-humped light curves. The long- 
term and short-term variability of this CV is discussed together with the spiral arm 
structure of an accretion disk in the context of observational effects taking place in 
bounce-back systems. 

Subject headings: stars: novae, cataclysmic variables — stars: dwarf novae — stars: 
brown dwarfs — stars: individual: SDSS J123813. 73-033933.0 



1. Introduction 



The object catalogued as SDSS J123813.73-033933.0 (hereaft er SDSS 1238) was i dentified with 
a faint (r = 17.82 mag) short-period cataclysmic variable (CV) by lSzkody et al.l (120031 ). The optical 
spectrum of SDSS 1238 shows a blue continuum with broad absorption features originating in the 
photosphere of a white dwarf surrounding double-peaked Balmer emission lines, formed in a high 
inclination accretion disk. T he orbital period of the system is P rb = 

0.05592(35)d = 1.34(l)/i, 

based on spectroscopic data (jZharikov et al.ll2006l ). The orbital period and the spectral features 
match those of WZ-Sge-type systems, but with this, similarities practically end. A number of 
observed aspects of the system differ from the majority of short period CVs. The most intriguing 
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characteristic, which we found in this system, is a sudden and fast rise in brightness up to ~ 0.45 mag 
during a short time, of about half of the orbital period. After reaching its peak, the brightness slowly 
decreases, lasting ~ 3 — 4 hours, down to the quiescence level. We call these events brightenings 
in order to distinguish them from the more common outbursts, humps, flickering and other types 
of variability documented in short period CVs. These brightenings seemed to happen cyclically 
about every 8-12 hours. In addition to brightenings, a nearly permanent sinusoidal variability was 
detected in the light curve of SDSS 1238 with a period half that of the spectroscopic orbital period 
P = -forb/2 = 40.25min (hereafter, the double-humped light curve). The amplitude of the double- 
hump variability depends on the phase of the brightenings. It increases with a total ri se in brightness 
of up to ~ 0.2 mag and decreases until almost disappears d u ring t he quiescence (jZharikov et al 



200 



y). A similar behavior was found later by Szkodv et al. ( 2006 ) in another short period CV, 



SDSS J080434.20+510349. 2 (hereafter SPSS 0804) wh ich has an identical spectral appearance to 
SDSS 1238 in quiescence. IZharikov et al.l (|2006l . 120081 ) advanced the hypothesis that the double- 
humped light curve is a signature of 2:1 resonance in the accretion disks of these systems. In order 
for the accretion disk to reach permanently the 2:1 resonance radius, the mass ratio of the binary 
component must be extreme {q < 0.1), and as such these objects could qualify as bounce-back 
systems, e.g. CVs, which are old eno ugh to rea c h the period minimum and leap toward slightly 



longer orbital periods, as predicted by iPaczynskil (|198ll ). It is supposed that accretion disks of WZ 



Sge systems reach 2:1 resonance radiu s during super-outbur st, when some of them have been noted 
to show double- humped light curves (jPatterson et al.ll2002j ). The super-outbursts of WZ Sge type 
systems are infrequent and happ en every two dozen or more years. Oddly enough, SDSS 0804 went 
into the super-outburst in 2006 (jPavlenko et al.l 120071 ) and exhibited all necessary attributes of a 
classical WZ Sge- type object. Regrettably, the brightenings disappeared from the light curves of 



SPSS 0804 after the super-outburst, although the double-hump light curve persists (jZharikov et al 



2008). Thus, the SPSS 1238 remains the only object that still shows brightenings. Intrigued by the 



new photometric phenomenon observed in these two systems, we conducted a new time-resolved 
photometric study of SPSS 1238 to establish the reasons behind their common nature, understand 
the origin of the cyclic brightenings and its relation to the amplitude of the double-humped light 
curve. Meanwhile, we discovered that SPSS 1238 was marginally detected as an infrared source by 
2MASS survey, and secured accurate near-IR photometry of the object. In Sect. 2 we describe our 
observations and data reductions. The data analyses and the results are presented in Sects. 3, 4, 5, 
while a general discussion and conclusions are given in Sect. 6. 



2. Observations and data reduction 

The object is listed in the 2MAsd3 Point Source Catalogue with J = 16.65(13), H = 16.49(23), 
K ~ 16.42 magnitudes. These magnitudes are close to the detection limits of 2MASS, particu- 



1 http://www.ipac. caltech.edu/2mass/ 
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larly in the i^-band. In order have more accurate photometry, we obtained new observations of 
SDSS J123813.73-033933.0 in J H K s on 17 June 2009 with the near-infrared camera PANIC (Mar- 
tini et al. 2004) attached to the 6.5 m Baade/Magellan Telescope at Las Campanas Observatory 
(LCO). PANIC provides an image scale of 0.125" pixel -1 on a Hawaii HgCdTe 1024 x 1024 array 
detector. The FWHM of the point-spread function was between 0.55" and 0.60" during our observa- 
tions. For each filter, 9 dithered frames spaced by 10" were taken, with total on-source integration 
times of 540 s in each of the three filters. The nine frames were shifted and averaged to produce the 
final images. Standard sky-subtraction and flat-field correction procedures were applied. Aperture 
(1.2") photometry was performed with DAOPHOT within IRAF in the stand ard way Flux calibra - 



tion was performed using standard stars SJ 9146 and SJ 9157 from the list of iPersson et al.l (|1998l ) 
and the total errors are estimated to be less than 0.05 magnitudes. The resulting magnitudes of 
SDSS 1238 are J = 17.07(5), H = 16.65(5), K = 16.42(5) and corresponding colors J -H = 0.42, 
H — K s = 0.23 and J — K s = 0.65. 

In order to investigate whether flux variations in the near-IR occur in this system in timescales 
of a few minutes, we measured the J, H and K s fluxes from each of the nine short-exposure (60 s) 
frames in each filter. In the ~ 15 min that lasted each series, we did not detect any variability within 
the ~ 0.15 mag photometric uncertainty associated with each single frame. Note, that in longer 
timescales (comparable to the orbital period), the system is expected to show some variability in 
the near-IR, mostly due to the elliptical shape of Roche-lobe filling secondary, but these could be 
missed in a 15 min time series. The ellipsoidal variability of the secondary can be calculated and 
it has been taken into account in further considerations of IR magnitudes. 

The objective of our optical photometry of SDSS 1238 was to study of phenomenon of bright- 
enings. Taking into account the long duration of the brightenings and the uncertainty of the cycle 



period (jZharikov et al.ll2006l ). we planned and executed a multi-longitude observational campaign 
of this object. Time-resolved CCD photometry was obtained at several facilities: the 1.5 m tele- 
scope at the Observatorio Astronomico Nacional at San Pedro Martir in Mexico; the 0.8 m IAC80 
telescope at the Observatorio del Teide in the Canary Islands, Spain; the 2.1 m telescope at the 
Bohyunsan Optical Astronomy Observatory (BOAO) in South Korea; the 2m telescope at the 
Terskol Astrophysical Observatory in the Northern Caucasus, Russia, the 2m Himalayan Chan- 
dra Telescope of the Indian Astronomical Observatory (IAO), Hanle, India. The data reduction 
was performed using both ESO-MIDAS and IRAF software. The images were bias-corrected and 
flat-fielded before aperture photometry was carried out. The log of photometric observations is 
presented in Table. 1. 

The long-slit observations have been obtained with the Boiler & Chivens spectrograph on 
the 2.1-m telescope at the SPM site with a resolution of 3.03 A pixel -1 . The spectra span the 
wavelength range 4000-7100 A. In order to improve the signal-to-noise ratio, we obtained a series 
of phase-locked spectra: 10 spectra were taken at equal phase intervals over a single orbital period 



2 www.astrossp.unam.mx 
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Porb = 80.5 min with an exposure time of 486 sec per spectrum. This sequence of spectra was 
repeated at exactly the same phase intervals for subsequent periods and subsequent nights. This 
allows us to calculate the phase- averaged spectra, summarizing the spectra of the same orbital 
phase obtained during one night and the whole set of observations without further decreasing the 
time resolution. The log of spectroscopic observations is presented in Table. 2. 



3. Spectral energy distribution: system parameters and distance to the object 



The spectral energy distribution (SED) of the object in the range of 4000-25000A is shown in 
FigJH The detailed d escription of the SDSS 1238 optical spectrum was given in our previous paper 
(jZharikov et al.l 120061 ). The overall appearance of the spectrum has not changed, but we detect a 
significant variability of the equivalent widths of Balmer emission lines from epoch to epoch. The 
Balmer lines are about two times weaker in the 2009 observations compared to the 2004 spectra 
(Figj2]). The continuum, however, has not changed during the last five years, as we compare V- 
band magnitudes in quiescence between brightening s. The average quiescence magnitud^l between 
brightenings remains constant at around V = 17.8 ± 0.1. 

The present near-IR measurements demonstrated that there is significant IR excess emission 
to that expected from the Rayleigh- Jeans tail of the optical spectrum for a white dwarf (WD). In 
fact, neither a WD nor a power-law flux from the accretion disk, nor their combination, can explain 
the observed IR excess. The most probable source of IR excess is the radiation from the secondary 
star. To determine the spectral type of the secondary and the distance to the system we fitted the 
observed optical-infrared spectral energy distribution of the object with a simple model: the total 
flux F*(X) is the sum of contributions from a W D with a hydrogen atmosphere, i*wD(Ieff> A) (DA 

7 J 1 1 I 

type WDs), an accretion disk with Fad ~ A~3 (jLynden-Belll 119691 ) . and a red/brown dwarf with 
Fbd(X): 

F*(X) = F WD (T eS , A) + Fad(A) + F|P T (A) (1) 

Brown dwarf fluxes were taken from the literature ( McLean et al. 2003 . 2007 ) and on-line sourcefl 
The white dwarf spectra with a mass range of M\vd = 0.6 — 1.1 M Q were used with a 0.1 M step 
and the radii were calculated using the white dwarf radius-mass relation of 



R WD = 1.12 x 10 y 1 



1.44M 



from lNauenbergl (|1972h and lWarnerl (119951 ). Spectra of WDs in th e 4000-25000A range with pure 
hydrogen atmosphere were obtained using ATLAS9 (jKurucall993l ) and SYNTH (|Piskunovi ll992) 



3 Secondary photometric standards were established in the field of SDSS 1238, by calibrating them using reference 
star S2003313360 from GSC-II www.gsss.stsci.edu/Catalogs/GCS/GSC2/GSC2.html 

4 see http://web.mit.edu/ajb/www/browndwarfs/ 
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codes for an appropriate range of temperatures. Although our previous temperatu re estimate was 



Twn = 15 600 ± 1000K based on fits to the absorption portion of Balmer lines (jZharikov et al 



2006), in the present fitting procedure we allowed a wider temperature range from T c g =11 000 to 
18 000 K, because of the larger number of free parameters. The calculations were performed with 
a 1 000 K step and with the surface gravity g = 7 p™ . The spectra are normalized to Ao = 5500A 

"WD 

and the contribution of the WD is 

F W d(T cS , A) = C x {5) * i^ m (T eff , A), 

where C\(5) = iq-° a *( v + s + m v) ; V = 17.8 is the object's magnitude in quiescence, and 5 is a 
parameter, determining the contribution of the flux from the WD in the V band. Finally, the 
My = 21.109 - is the constant to convert magnitudes into flux (in ergs/cm 2 /s/A) in the V band. 
The spectra of the accretion disk was assumed to be a simple power law 

7 

Fad (A) = (Ci(0)-Ci(<5)) x (A 

where (Ci(0) — Ci(5)) determines the contribution from the accretion disk in the V band, assuming 
that the WD and the accretion disk are the only contributors in that wavelength as the only other 
contributor is the brown dwarf, which has a negligible flux in V. 

The distance to the object is estimated to be 

d = i?WD 



F bb (T eff ,5500A) 



F WD (5500A) ' 

where F hh (T e g ; 5500A) is the black body flux at A = 5500A with effective temperature T e ff . 

Observed SEDs of red/brown dwarfs of spectral types between M6 to L5 normalized to fit the 
observed flux i n J were used. The b olometric correction to the J magnitude for each spectral type 



was taken from lTinnev et al .1 (|2003l ). 



The free parameters of the three-component model are: the white dwarf effective temperature, 
T e fj , the mass of the white dwarf, JV%d, the spectral type of the secondary star, SpT, and the 
parameter 5 (in magnitudes). The best fit model to the observed SDSS1238 spectrum in the 0.4- 
2.5/um range achieved for the following set of parameters (FigfT]) is: Twd = 12000K, SpT = L4, 
M^D = 1.0M Q and 5 = 0.1. The deduced distance to the object is 110 pc. We studied the behavior 
of x 2 vs a single fit parameter, when the other three are fixed to their corresponding best values. 
Figj3] presents x 2 plots for various parameters with marked confidence levels corresponding to 95, 
80 and 60%. We can state that x 2 tends to the minimum value always, when Twd = 12 000 K 
and 5 = 0.05 — 0.15 regardless of the value of the other parameters. At the same time, there is 
a dependence between the mass of the WD and the spectral type of the secondary: the lower the 
mass, the larger the radius of WD, and thus, the larger the distance to the system, resulting in an 
earlier spectral type of the secondary. The best fit to the optical part of the spectrum is reached 
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with the mass of the WD of Mwd = 1-0 Mq, leading to the cited distance of, 110 pc, and spectral 
type of the secondary, L4. The entire range of secondary from SpT = M9 at Mwd = O.6M and 
d = 160pc to SpT = LA at Mwd = 1.0 Mq, d = HOpc was considered. However, a pronounced 
minimum in \ 2 for rather massive WD strongly suggests the presence of brown dwarf in this system. 
The indication of a massive WD is not accidental, a CV with a brown dwarf secondary, supposed 
to have evolved beyond minimum orbital period limit, has an age ~ 3 — 5 x 10 9 years and a long 
history of accretion. A priori, the bounce-back CVs are assumed to harbor a massive WD, a lthough 
the d istribution of known masses of all WDs in eclipsing systems does not s how any trend ([Kn igge 
2006). But when considering only the short-period end of that distribution as lLittlefair et al.l (|2008) 
did, then it becomes apparent that evolved systems systematically have more massive WDs. The 
upper limit for the mass of th e secondary is Mbd < 0.09M Q for SpT = M9 or Mbd < 0.07M Q for 
SpT = LA (jClose et al.l 12003). But it is also a well known fact that the secondaries of CVs tend 
to show systematically an earlier spectral class and a lar ger radius with secon daries of lower mass 
than the corresponding single stars on the main sequence (jLittlefair et al.ll2008l ). So we expect that 
the mass of the secondary in SDSS 1238 hovers around the lower edge of the above mentioned range 
of masses. 

Note, that due the high galactic latitude of this object (b = 59.5) and the inferred small 
distance, the interstellar extinction is negligible. We also would like to emphasize that in the 
above calculations we took into account the fact that the state of the system was unknown at the 
moment of acquiring the IR data. Therefore, we conducted the fitting for both cases, considering 
that the system might have been at the maximum of the brightenings during the IR observations 
or at the bottom. That introduced only a minor change, basically decreasing the distance to the 
system by ~ 15pc and not affecting our conclusions regarding the spectral type of the secondary. 
Finally, we would like to comment on the discrepancy in th e temperature determ ination of the 
WD, which is larger than we would like as compared to that of Zharikov et al. (j2006l ). In the latter 
the temperature and gravity could not be determined simultaneously, and one parameter had to 
be fixed in order to calculate the other, which always introduces ambiguity, as none of these two 
parameters could be estimated independently. In this paper we take a more complex approach: 
not only profiles of the lines are being fitted, but the continuum is taken into consideration as well. 
Also, the presence of the secondary, adding additional restraints on the distance and thus, on the 
size of the WD. The SED of the accretion disk, which may not necessarily obey the canonical power 
law index is the source of the largest uncertainty in our analysis. Since the contribution of the disk 
is limited to only < 20 %, so is the accuracy of our estimates. With all that in mind, we still end 
up with a range of spectral classes for the secondary that implies that SDSS 1238 is a bounce-back 
system. 
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4. Doppler tomography of SDSS1238 



The Hq, emission line originating in the accretion disk is the least affe cted by the absorption 



from the underlying white dwarf. Therefore, we constructed Doppler maps (jMarsh fc Hornd ll988) 
of H Q using all our data obtained in 2004 and 2009 (FigJU top). The individual maps of separate 
nights resemble each other so much, that combining all available data did not smooth out details, 
but made them much more convincing. In order to over-plot contours of the secondary star, the 
location of the WD, the trajectory of the stream and resonance radius of the disk, we used the best 
fit parameters to the SED (see previous section). Thus, the white dwarf mass of Mwd = 1-0 M & was 
adopted. Due to the ambiguity in determining precise masses of brown dwarfs, we adopted a mass 
ratio of q = 0.05. This value is typical for systems considered as bounce-back (jKniggel 12006) and 
it is the same mass ratio which was obtained for SDSS0804, a twin of SPSS 1238, from the super 



hump period observed during the super-outburst in 2006 (jPavlenko et al.l 120071 ; IZharikov et al 



2008). Before SDSS0804 underwent a WZ-Sge type super-outburst, it showed similar peculiar 
photometric variability to SDSS 1238. In addition, it shares every other charact eristic of a WZ 



Sge-t ype object, and was also proposed as a candidate to the bounce back system (IZharikov et al 
200sh . 



As already noted, the structure of the accretion disk did not change between the two epochs 
of our observations. There is a bright spot at the expected place where the stream of matter from 
the secondary collides with the accretion disk, but it overlaps with a much larger and prolonged 
structure, too extended to be a part of the spot. Another extended bright region of similar size can 
be seen at velocity coordinates (~ 700 km s , « km s -1 ) as well as a less bright structure at 
(~ —200 km s — 1 , ~ —800 km s" 1 ]. Similar Doppler maps were obtai ned for WZ Sge during super- 



outburst in 2001 ( Baba et al 



2002; 



Howell et al 



back candidates SDSS 1035 (jSouthworth et al.ll2006l ) and SDSS 0804 (lAviles et alJIin preparation! ). 



2003; ISteeghsl I2004J ) and in quiescence for bounc e 



Such brightness distribut ion in the Doppler map c an b e interpreted as evidence of spiral waves 
in the disk (see for example, ISteeghs fe Stehld (|1999l ) and ISteegha (|200ll ) and reference therein). 
Th e formation of a spiral s tructure in an accretion disk of a close binary sy s tem w a s predicted 



1979) and explo r ed by various auth ors ( Matsuda et al. ( 199oh l 



Heemskerk 



Kunze &: Speithl (J2005J) ; iTrussI (120071 ) and reference therein) . ISawada et al 



bv iLin Sz Papa oizou 
( 1994h:|stehlel jl99flh : 

(|l986al lbh demonstrated from high resolution numerical calculations that spirals will always formed 
in accretion disks under tidal forces from the secondary. They actually used q = 1 in their models, 
but observationally, such spirals were detected in a number of systems only during outbursts of 
dwarf novae. The careful examination of quiescent disks of th e same systems did not reveal any 
spiral structures in longer period DNe. ISteeghs &; Stehld (|1999l ) argued that little evidence of spiral 
arms in the emission lines is expected in systems with low values of viscosity. 

On the other hand, spiral arms related to 2:1 reso nance can be found i n syst ems with extremely 
low mass ratio q < 0.1 as originally was predicted by lLin &; Papaloizoul (|1979j). The bounce-back 
systems and related to them, WZ Sge stars, are examples of such objects. The long outburst recur- 
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rence time in WZ Sge systems is probably explained by a very low viscosity in their accretion disks, 
yet spiral arms can be observed permanently in quiescent bounce-back systems in which, on one side 
there is a massive WD, which gained mass during a long accretion history, and on another side there 
is a late-type brown dwarf, providing a mass ratio of < 0.06. Figjl] (bottom, left) depicts a synthetic 
Doppler map constructed from a model accretion disk that is shown on the bottom right panel. The 
latter was calculated with t he binary system parameters derived above by using smoo th particle 
hydrodynamics according to iMurrayi (|1996l ) ; iKunze et al.1 (|1997. and references thereinl ) . The arti- 
ficial doppler map reproduces the observed map in a case when there is a brightness excess within 
spiral arms. Most of the disk particles are on periodic orbits, which are most favorable from the 
point-of-view of viscosity. However, the resonance dispatches some particles onto aperiodic orbits 
creating viscosity perturbations, whic h will create ex cess of heat. A slightly different interpretation 
of spiral arm brightness is offered by lOgilvid (|2002l ). The mechanism is not very well established, 
but it is natural to assume that in these regions there will be excess emission. 



5. Cyclic brightenings. 

Figure [5] displays the light curve of SDSS 1238 ob tained in 2007. I n gen eral, the object shows an 
identical behavior to previous years, as described by Zharikov et al. (j20od ). There are two distinct 



types of variability: a long-term variability (LTV), lasting more than 8h, and a short-term variability 
with a period corresponding to half the orbital period. Follow-up observations during 2008 and 
2009 confirm a steady presence of both types of variability in the light curve (Figj6]). Continuous 
observations during about 15 hours obtained on 22/03/2007 (HJD 2454182.15-2544182.60, FigJT]) 
allowed us to obs erve two consecutive brightenings and, thus to determine their recurrence time 



directly (~ 9h). I Zharikov et al.l (|2006l ) demonstrated that brightenings occurred cyclically with 
periods in between 8 and 12 hours, however it was not possible to establish if the phenomenon 
was strictly periodic or not. Armed with more data, and the advantage of detecting consecutive 
events by employing multi-longitude observations, we performed a simple period search. T he period 



analy sis of all data obtained during 2007, based on a discrete Fourier transform method (IDeeming 



19751 ). results in a strong peak at a frequency /ltv = 2.59 day -1 with FWHM of 0.1. This conforms 
with the period of Pltv = 9.28±0.36h « 7P or b (FigJHJ). However, the data folded with this period 
looks messy because, apparently, some re-brightenings happen with a different cycle period. It 
shows that the large uncertainty in the period value is not just a result of a scarce amount of data, 
or its uneven distribution, but that the brightenings are not strictly periodic in long timescales. 
The peak 2 x / or b = 35.76 day -1 , corresponding to the half of the orbital period, -P rb/2 = 40.3 min, 
is also present in the power spectrum. There are some additional peaks which are one-day aliases 
or a combination of high harmonics of /ltv with 2 x / or b- We have repeated the same analysis 
for the data obtained only within HJD 2454177.1-2544185.6, when we observed the object with 
small observational gaps. The result of the analysis of periodicity for this selection is presented 
in FigEJ The strongest peak in the power spectrum appears at a frequency corresponding to 
P£ TV = 9.34(26)h. The light curve folded with this period is coherent during the considered time 
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(Figj9]), i.e macro-profile of the brightenings (the rise, the peak, the declining slope and minimum 
duration) repeats itself with high accuracy and, if not for the high-frequency modulation also 
present in the light curve, the folded curves would concur. Combining the entire available data 
obtained during the 2004 — 2009 observational runs does not permit a determination of a single 
period with which the entire data set can be folded. Clearly, the LTV modulation is not strictly 
periodic, but probably maintains coherence in a time scale of several weeks and certainly has a 
cyclic nature. 

The long-term variability (the brightenings) presents a fast rise in the object brightness of up 
to 0.4-0.6 mag, lasting about half the orbital period (~ 0.7 h) with the brightness falling back to 
the quiescence level during the next 4.4-4.6 h (Fig. [7]). After that, the object remains in quiescence 
during about another 4 hours. The brightness increases at a rate of ~ 0.75 mag/h, and the falling 
rate is about ~ 0.10 mag/h. Sometimes, the rise in the brightness happens more slowly and lasts 
until almost one full orbital period (Fig[5] nights: 156, 179). 

A short-term variability is clearly present during the brightenings. The amplitude of the short- 
term variability depends strongly on the phase of brightenings. The amplitude is larger then the 
total brightness increases and falls, sometimes practically disappearing, between the brightenings 
(Figs. U II Ef. Curiously, when the LTV is coherent, the signal at double orbital frequency 
decreases and a beat frequency between the LTV and 2 x / or b becomes dominant in the range 
of frequencies around 36 cycles/day. Removing LTV has little effect, because this does affect the 
frequency modulation but not the amplitude modulation. The strong relation between LTV cycle 
coherence and the strength of 2 x / or b modulation attests that these two phenomena are physically 
linked in more than a simple amplitude correlation. We suspect that the brightenings and formation 
of spiral arms in the disk are the result of the same process, but its cyclic nature remains unclear. 



6. Discussion and conclusions 

In recent years, a number of objects were discovered, mostly thanks to the Sloan Digital Sky 
Survey, which can be characterized as WZ Sge type CVs in quiescence, judging by their ~ 80 — 85 
min periods, spectral appearance and lack of SU UMa style outburst activity. Currently, there are 
about 20 objects with similar characteristics which were selected from the SDSS lists. Three such 
objects are GW Lib, V455 And and SDSS0804, which underwent a typical WZ Sge type super- 
outburst in 2006-2007, confirming the correct assessment of their nature from their behavior in 
quiescence. But unlike classical WZ Sge stars, some of them differ by the presence of permanent 
double-humps per orbital period in the light curve. The significance of the double-humped light 
curve lies in the assumption that it is produced by the spiral structure in the accretion disk. The 
spiral arms in accretion disks have been detected before by means of Doppler tomography in long 
period systems only during DN outburst. We are not aware of any reports of double-humped light 
curves during these detections. Double-humped light curves have previously came to attention, as 
they were detected during super-outbursts of WZ Sge. It was soon suggested that they are result 
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of 2:1 resonance (jOsaki fc Meverl 120021 ; iPatterson et al.ll2002l ). as one of the possibilitie s. Also 



spira ls arms are formed in the accretion disk of the WZ Sge undergoing super-outburst (jSteeghs 
2004 ). because, according to our hypothesis WZ Sge accretion disk reaches the resonance radius 
only during super-outburst. Nevertheless, the spiral structure has been observed persistently along 
with the double-humped light curve in two similar short-period systems, SDSS 0804 and SDSS 1238 
in quiescence. We assume that the spiral structure formed in low viscosity, quiescent disks is a 
result of distortion of the disk by the 2:1 resonance. The 2:1 resonance may happen only in the 
accretion disk of a system with extreme mass ratio of q < 0.1. Such mass ratio is achieved only 
in systems known as bounce-backed or, in other words, systems which hav e reached a minim um 
period of sa 80 min and have turned to slightly lon ger periods according to iPaczvnskil (|198ll ) . It 
is expected that bounce-back systems are numerous (jKolb Baraffdll999l ). but until recently very 
few candidates have been found. The SPS S helped to uncover a large number of new CVs, with new 
interesting features (jGansicke et a l. 2009). Among them, there is a number of short period systems , 



some of which turned out to belong to the long-sought bounce-back systems (jLittlefair et al 
Mennickent fc Diazll2002h . 



2008 



We have unveiled a brown dwarf secondary in SDSS 1238, probably as late as L4, thus providing 
strong evidence that this object is a real bounce-back system. We also estimated parameters of the 
WD. The best fit to the SED converges if the primary is a massive Mwd ~ 1-0 Mq white dwarf 
with a T\vd = 12 000 K temperature. Both numbers seem plausible since the system is very old 
and the WD is expected to be relatively massive and cool. These findings provide further support 
to our claim that by simply observing a permanent double humped light curves one can identify 
an evolved, bounce-back systems instead of using other complicated methods. 



It is important to note that SDSS 1238 has another peculiarity, which IZharikov et al.l ([2006) 
termed as brightenings. There was only one other system known to exhibit brightenings, namely 
SDSS 0804, but since it underwent a super-outb urst in 2006, i t s phot ometric behavior has drastically 
changed. The brightenings shortly detected bv lSzkodv et al.l (|2006l ) in SDSS 0804 before the super- 
outburst have been replaced by: a) the mini- outb urst activity with pe rmanent presence of the 
double-humped light curve of constant amplitude (IZharikov et al.1 I2008T) and by: b) a 12.6 min 
period, probably corresponding to pulsation activity of the WD (jPavlenko et al.1 120071 ). Therefore, 
at present, SDSS 1238 is the only object known to show brightenings. It greatly complicates the 
study of this phenomenon. Based on new multi-longitude continuous monitoring, we demonstrated 
here that the brightenings are of cyclic nature with a recurrence time of ~ 9 hours and they are 
probably coherent over several cycles. There seem to be strong correlation between brightenings 
cycles and the amplitude of double-hump periodic variability. This is believed to be the result of 
spiral arms in the accretion disk of bounce-back systems with an extreme mass ratio, in which 
accretion disk extends beyond a 2:1 resonance radius. The tomogram of a simulated accretion disk 
in the regime of resonance, closely resembles the observed one and supports this hypothesis. The 
formation of a spiral structure in the disk can be accounted by the appearance of double humps in 
the light curve, but it can not be the reason of increased brightness of the disk. The disk brightness 
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directly depends on the mass transfer rate and its change should probably reflect change in mass 
transfer. No readily explanation is available as to why that rate can be variable and cyclical, but 
possible speculations include counteraction: a) to the heating of the secondary by brightenings, or 
b) to tidal interaction between the secondary with the resonance attaining accretion disk. 

This work was supported in part be DGAPA/PAPIIT projects IN109209 and IN102607. 
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Table 1: Log of time-resolved observations of SPSS J123813.73-033933.0 in V band 
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Table 2: Log of time-resolved spectroscopic observations of SDSS J123813. 73-033933.0 



Date 
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Fig. 1. — The spectral energy distribution of SDSS 1238 and the result of the model fit (top). 
Spectrum fragments around H 7 , EL; and H Q lines are shown in bottom panels. 



-18- 




Fig. 2. — The low-resolution time-average spectra of SDSS1238 obtained in different epochs. 
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Fig. 3. — x 2 vs parameters of the fit, where T e ff is an effective temperature of the primary 
WD, M\vd is a mass of the WD, SpT is a spectral type of the secondary (from M6V/SpT=6 to 
L6/SpT=16) and 5 is a ratio between accretion disk and WD contribution in the continuum flux of 
the SDSS1238 in the V band. The formal confidence levels by the fit are presented by dashed lines. 
The numbers at the dashed lines are a probability to reject a model with \ 2 above corresponding 
line. 
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Fig. 4. — top) H Q Doppler maps constructed on all data obtained in the 2004 (left) and 2009 (right) 
yy. The circle show the velocity at 2:1 resonance radius, bottom) The synthetic Doppler map (left) 
obtained from a model accretion disk (right) for the system with yiwD = 1-0Mq and q = 0.05. 
The circles correspond to 2:1 and 3:1 resonance radiuses. 
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Fig. 5. — The light curve of SDSS 1238 in V band obtained in the 2007 year. Each night is presented 
in a separate panel. 
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Fig. 6.— The light curve of SDSS 1238 in V band obtained in 2008 and 2009 yy. Each night is 
presented in a separate panel. 
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Fig. 7. — The light curve of SDSS 1238 in V band obtained during about 15 hours continuos obser- 
vations in HJD 2454182. Vertical dashed lines select different phases of the long-term variability 
(LTV). The numbers in the top of the figure are approximately durations for each marked phase 
of LTV and the numbers in the bottom are corresponding rates of the magnitude change during 
fading and increasing of the object brightness. 
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Fig. 8. — top) The power spectrum of the all photometric data presented in the top panel of Figj6j 
bottom) The light curve comprised of the all data folded on P = 9.28h. The open color points 
show brightenings with maximal displacment in the light curve folded on P = 9.28h. 




Fig. 9. — top) The power spectrum of the data obtained in the period of HJD 2454177.1-2454185.6 
(see FigE]). bottom) The light curve comprised of the selected data from the period of HJD 
2454177.1-2454185.6 folded on P = 9.34h period. 



